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ABSTRACT 

The short GRB 120323 A had the highest flux ever detected with the Gamma-ray Burst Monitor 
(GBM) on board the Fermi Gamma-ray Space Telescope. Here we study its remarkable spectral 
properties and their evolution using two spectral models: (i) a single emission component scenario, 
where the spectrum is modeled by the empirical Band function (a broken power law), and (ii) a 
two component scenario, where thermal (a Planck-like function) emission is observed simultaneously 
with a non-thermal component (a Band function). We find that the latter model fits the integrated 
burst spectrum significantly better than the former, and that their respective spectral parameters are 
dramatically different: when fit with a Band function only, the Ep^ of the event is unusually soft for 
a short GRB (70 keV compared to an average of 300 keV), while adding a thermal component leads 
to more typical short GRB values (Epeak ^300 keV). Our time-resolved spectral analysis produces 
similar results. We argue here that the two-component model is the preferred interpretation for 
GRB 120323A , based on: (i) the values and evolution of the Band function parameters of the 
two component scenario, which are more typical for a short GRB, and (ii) the appearance in the 
data of a significant hardness-intensity correlation, commonly found in GRBs, when we employee 
two-component model fits; the correlation is non-existent in the Band-only fits. GRB 110721A, a 
long burst with an intense photospheric emission, exhibits the exact same behavior. We conclude 
that GRB 120323A has a strong photospheric emission contribution, first time observed in a short 
GRB. Magnetic dissipation models are difficult to reconcile with these results, which instead favor 
photospheric thermal emission and fast cooling synchrotron radiation from internal shocks. Finally, 
we derive a possibly universal hardness-luminosity relation in the source frame using a larger set of 
GRBs (LP and = (1.59 ± 0.84). 10 50 (E^ )i ) 1 - 33±0 - 07 erg.s -x ), which could be used as a possible redshift 
estimator for cosmology. 

Subject headings: Gamma-ray burst: individual: GRB 120323A - Radiation mechanisms: thermal - 
Radiation mechanisms: non-thermal - Acceleration of particles 



1. INTRODUCTION 



Gamma-Ray Bursts (GRBs) are extremely energetic 
explosions at cosmologica l distances ( Meegan et al-11992 1 



van Paradijs et al. 1997) resulting most likely from the 



formation ot stellar mass black holes, either through the 
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collapse of massive stars (Wooslcy 1993; MacFadyen & 
Woosley|1999 Woosley &: H cger 2006 1 or v i a the merger 



ot two compactTobj ects QPaczynskT 



1986 Fryer et al. 



1999 Rosswog 2003). Regardless of the nature and for 



mat ion mechanism of the GRB central engine, the fire - 
ball model, initially proposed by Cavallo & Rees (1978), 
best explains the GRB emission from radio up to Gev 
gamma-rays. During the GRB explosion, the central en- 
gine produces a collimated bipolar jet, mainly composed 
of electrons, positrons, photons and a small amount of 
baryons. The central engine ejecta are accelerated to 
rclativistic velocities forming layers of high density re- 
gions, which propagate at various speeds. When the 
faster layers catch up with the slowest, charged particles 
are accelerated through mildly relativistic collisionless 
shocks: this is the so-called internal shock phase (Rccs 
fc Meszarosl | 1994| |Kobayashi et al.| [1997{ |Daigne & 



Mochkovitch 



1998). As a result they then produce non 



thermal radiation such as synchrotron emission, observed 



(see e.g., the si 


jectral catalogs by Preece et al. 


2000 


Goldstein et al. 


2012a), and even up to several tens o 


GeVs (see e.g., 


Abdo et al. 2009a |Abdo et al. 


2009b 


Ackermann et al. 2010 2011). As the jet expands, it 



in a relativistic shock; the charged particles involved in 
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this collision produce synchrotron radiation visible from 
radio wavelengt hs to X-ray energies . This is the exter- 



nal sh ock phase (Rees & Meszaros 1992 Meszaros & Rees 
|1993| , which is responsible tor the afterglow emission ob 



served during few hours to several days and even years 
following the prompt emission. Alternative models for 
the GRB prompt phase are magnetically driven, involv- 
ing mechanisms such as magnetic field line reconnection. 

Besides the non-thermal (synchrotron) radiation, the 
fireball model also predicts strong thermal emission em- 
anating from the jet's photosphere, which would be ob- 



Thompson scattering ( 


Goodman 


1986 


Meszaros 


2002 


Rees & Meszaros |2005 


). In the most standard version 



erated outflow (fireball), this thermal emission would be 
very int ense and would overpower the n on-thermal com- 
ponent (IDaigne fe Mochkovitch||2002[ ). Prompt GRB 
spectra nave been traditional ly adequately fitt ed with 
the empirical Band function (Band et al. 1993), which 
is a smoothly (with curvature) broken power law, with 
indices a and (3 for the low and the high energy part, 
respectively. The break energy of the Band function, 
parameterized as Ep^ i , correspond s to the maximum 
of the v¥ v spectrum (Gehrels 1997), when a >-2 and 
[3 <-2. The Band parameters are usually compatible 
with non-thermal emission. The synchrotron mechanism 
thus remains the preferred model to explain most of the 
prompt emission; note, however, that the a parameter 
values are often incompatible with slow and fast electron 



cooling scenarios ( |Crider et al.]|1997| [Preece et al.||1998[ ). 

However, a few bursts observed with the Burst 



And Transient Source Experiment (BATSE) onboard 
the Compton Gamma Ray Observatory (CGRO) were 
found to be well fitted by a single blackbody func- 
tion, th e tell-tale signature of emission f rom the photo- 
sphere |^irlanda^erarj2003j [Ryxte]2004). Furthermore, 
time resolved analysis of the strongest GRB pulses ob- 
served with BATSE were shown to, in some cases, be 
well fitted by a combinat ion of a blackbody and a power 
law function ( |Ryde||2005[ ) . The pulse temperatures were 
found to lie within 20 — 100 keV and were observed to 
evolve in a c haracteristic way, dec aying as a broken power 
law in time ( |Ryde fc Peer|[2009| . However, the BATSE 
energy range (25 — 1800 ke V) made it difficult to fully 
assess these models (see e.g., |Ghirlanda et al.|2 007). The 
broader energy range (8 keV - 40 MeV) of trie Gamma- 
ray Burst Detector (GBM) on the Fermi gamma ray 
space telescope a l leviat es this shortcoming. Using Fermi 
data, Ryde et al. (2010) proposed that the dominant con- 
tribution to the spectru m of GRB 090902B is a modified 
blackbody. Furthermore, Guiriec et al. (2011a I reported 
for the first time the simultaneous fit to the data of long 
GRB 100724B using a blackbody (BB) and a Band func- 
tion corresponding to the thermal and non-thermal emis- 
sions, respectively. In this GRB, the thermal emission 
identified in the time-integrated spectrum and followed 
in the time-resolved spectroscopy analysis was a subdom- 
inant contribution, corresponding to only few percent of 
the total energy. Interpreted as a photospheric emission, 
this low intensity BB did not support the standard fire- 
ball scenario, where the energy initially released by the 
central engine is only thermal; instead, it suggested an 
initially magnetically dominated outflow (note that this 



conclusion is independent of the nature of the mechanism 
responsible for the non-thermal emission, internal shocks 
or magnetic reconnection). 

Another example is given by the strong emission pulse 
in GRB 110721A in which the presen ce of two com- 



ponen ts is yet again highly significant (Axelsson et al. 
2012|. The temperature of the BB component is ob- 



served to have the same characteristic temporal evolu- 
tion as seen in some of the BATSE bursts (Rydc 2005). 
On the other hand, the peak energy of the Band compo- 
nent decreases more rapidly and as a single power law in 
time. In the case of GRB 110721A the flux contribution 
of the BB is at most 10 %. 



See also Bhat & Guiriec 
review of recent GRB resu. 



ts. 



(2011) for an comprehensive 



Here we discuss the possible identification of such a 
subdominant but intense BB component in the short 
GRB 120323A. In section [2] and [3] we describe the Fermi 
observations of GRB 120323A and the analysis proce- 
dure, respectively. We present the results of our time- 
integrated and time-resolved spectral analyses in sec- 
tion H and [5l respectively. In section [6j we examine the 
hypothesis that the BB component might be present to- 
gether with a Band component during the entire burst 
and compare these results against Band fits only. In 
Section [7] we present an intriguing result on the relation 
between E p ° c s ^ k ; and the Band function luminosity, which 
supports and extends the often discussed E pca k-L correla- 
tion as well as the scenario of the simultaneous existence 
of the thermal and non-thermal components. We present 
a discussion and interpretation of our results in section [S] 
and conclude in Section |9| Appendix [A] describes the 
simulation procedure usedto assess and validate our re- 
sults. 

2. OBSERVATIONS 

The GBM onboard Fermi detected a very intense short 
burs t, GRB 120323A, on 2012 Ma rch 23 at 12:10:19.72 
UT ( |Gruber fc Connaughton1[2012| ). GRB 120323A has 
the highest peak flux among all events observed with 
GBM thus far; the intensity of the burst fulfilled the cri- 
terion for an Autonomous Repointing Request (ARR) of 
the Fermi spacecraft to place the source in the field of 
view of the Large Area Telescope (LAT). However, the 
burst was unusually soft and was not det ected at high en- 



ergies by the LAT ( |Tam fc Kongp012[ ) in the standard 
LAT data (100 MeV to >300 GeV), nor in the Low LAT 
Energy (LLE) data (designed to increase the LAT event 
acceptance at low energies and enable spectral analysis 
below 100 MeV). Figure [l] shows the GBM light curves 
of GRB 120323A in four energy ba nds ranging from 8 
keV to 1 MeV. The T90 duration (|Kouveliotou et al. 



1993|) of the e vent was T 90 = 0.448 ±0 090 s between 50 
an7T300 keV ( |Gruber fc Connaughton|2012| ). The evolu- 
tion of the spectral lag as a function of energy is shown 
in Figure [2j The spectral lags are calculated over the 
full duration of the GRB (T -0.02s to T +0.68s) using 
the cr oss-correlation method described in INorris et al.l 
(|2000|). The maximum spectral lag fo r GRB 120323A is 



small as expected for s hort bursts (INorris & Bonnell 
2006 Zhang et al. 2006). While the energy dependence 



of the lags is also consi stent wit h that for other intense 
short bursts ( Guiriec et al.|2 010), it is important to note 
that this may be due to pulse confusion between the two 
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Fig. 1. — Background substracted count light curves of GRB 120323A in four energy bands ranging from 8 keV to 1 MeV as observed 
with GBM. The light curves from 8 to 300 keV are obtained after summing the signal detected in Nal detectors NO, Nl and N3, and the 
light curve from 300 keV to 1 MeV corresponds to the BGO detector B0. The figures on the left show t he 2 ms light curves while the right 
ones correspond to the light curves rebinned for the fine time-resolved spectroscopy describe in section [5. 2| 
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Fig. 2. — Spectral lags measured between the lowest energy band 
8-20 keV and higher energy bands in the Nal and BGO detectors 
for GRB 120323A plotted as a function of mean energy of the 
higher energy band. The higher energy bands are: for Nal 20-40 
keV, 40-70 keV, 70-140 keV, 140-270 keV, 270-525 keV, 525-985 
keV; for BGO, 110-250 keV, 250-550 keV and 550-1380 keV. 



main emission peaks, with the second peak being spec- 
trally harder than the first and so contributing more to 
the cross-correlation function peak at higher energies. 

GRB 120323A was also detected with Konus- 
WIND ( |Golenetskii eFaT1[20T2] ) and MESSENGER. The 
best location tor this event was estimated with the Inter- 
Planetary Network (IPN) using all three satellites, to be 
centered at RA=340.4°and Dec=29.7°, inside an irreg- 
ular error box with a maximal dimension of 0.75° and a 



minimal dimension of 0.25° (Golenetskii et al. 2012 1. 



3. DATA ANALYSIS PROCEDURE 

We performed a spectral analysis of GRB 120323A us- 
ing only GBM Time Tag Event (TTE) data (8 keV - 40 
MeV). In the TTE data, each event detected with a GBM 
detector is recorded with its trigger time and the detec- 
tor energy channel. TTE data have the finest time and 
energy resolution and are then ideal to perform spectral 
analysis of such short GRB. For more information about 
GBM, see IMeegan et al.| p009|. We selected the three 



Nal detectors, 1N0, INI and IN 3, with angles to the source 
below 50° as these are not affected by blockage from other 
parts of the spacecraft nor shadowed by other detectors. 
We also used one BGO detector, B0, with a direct view 
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to the source. Since there is no detection of this GRB in 
the regular or LLE LAT data, we did not include these 
datasets in our spectral analysis; however, we note here 
that the extrapolation of our GBM spectral analysis is 
consistent with the LAT upper limits. 

We selected the Nal energy channels from 8 keV to the 
overflow channels starting at ^900 keV, and the BGO 
data from 200 keV to the overflow channels starting at 
~40 MeV. We then generated the response matrices for 
each detector using the best known location for the event, 
which is the IPN location reported in Section [2j For each 
detector, the background was estimated by fitting a poly- 
nomial function to time intervals pre and post burst. The 
background during the GRB was then estimated by ex- 
trapolating the polynomial function over the source time 
interval. Data were fit using the spectral analysis pack- 
age Rmfit provided by the GBM instrument team; the 
spectral fit was performed using the full spectral resolu- 
tion of the instruments. We determined the best spectral 
parameters by optimizing the Castor C statistic value. 
Castor Cstat (henceforth Cstat) is a likelihood technique 
modified for a particular data set to converge to a \ 2 with 
an increase of the signal. 

4. TIME-INTEGRATED SPECTRAL ANALYSIS 

We first analyzed the spectrum of GRB 120323 A over 
the entire duration of the GBM prompt emission, from 
T -0.016s to T +0.548s. In Tableland in Figure [3] we 
report the results for the various acceptable combinations 
of the models we tested. We discuss our results below, 
progressing from single to multiple component fits. 

A single component Band function fit gives a value 
for Ep Ca k of ^70 keV, which is in the tail (2-3%) of the 
Epeak distribution for GRBs observed with GBM and 



BATSE ( |Goldstein et al.||2012a| |Paciesas et~aT||1999| >. 

This low value ot E pea k is very unusual tor intense short 
GRBs, whose E pea k values are typically much higher than 
those of long ones (Paciesas et al. 1999 Guiriec et al. 
20101. A high redshilt (z) tor this GRB could reconcile 
the observed low E poa k value with the typical E pea k dis- 
tribution. However, while Epoak evolves as (1+z), the 



luminosity evolves as 47rDj 



Thus, a large distance 



would also increase dramatically the intrinsic luminos- 
ity of this GRB which is already extremely high. The 
energy flux (8 keV - 1 MeV) of GRB 120323A in the 
observer frame computed from Tq-0.016s to Tq+0.548s 



is (1.95±0.02)xl0 -5 e rg cm^.s 



which is above all 
In ad- 



the values reported in Goldstein_et al. (2012a 
as we show in Section [7J a redshitt of 
/itn 



'1 for 



dition, 

GRB 120323A is compatible with the observations. 

The Band function alone is not sufficient to describe 
adequately the time-integrated spectrum. We find that 
more complex models significantly improve the fits. A 
double broken power law (2BPL) and a Band function 
with a cutoff in the high energy power law (B+Cut) give 
the largest Cstat improvement over Band alone (60 and 
49 units, respectively, for two additional degrees of free- 
dom - dof ) , which suggests the existence of two spectral 
breaks, one around few tens of keV and the other around 
1 MeV. 

We next used a combination of a Band function with 



a BB (B+BB) as proposed in Guiriec et al. (|2011a) and 




Fig. 3. — Time-integrated spectral analysis from To-0.016s to 
To+0.548s. The figures on the left side show count spectra and 
the figures on the right show the corresponding model in the uF u 
space. The models used are from top to bottom, Band, B+Cutoff, 
2BPL, and B+BB. For each figure, the lower panel correspond 
to the residuals of the fit. For presentation purposes we rebined 
the energy channels in the figures, although the fit procedure is 
performed using the full spectral resolution of the instruments. 



BB affects the parameters of th e Band function simil arly 
to what was already reported in Guiriec et al. (2011a) for 
GRB 100724B: both a and j3 are shifted towards lower 
values, and E poa k changes from ^70 keV to ^200-300 
keV, a more typical value for a short GRB. We also notice 
that the temperature of the BB, kT ^10—13 keV, cor- 
responds to a vF v spectrum with a maximurr| 12 | around 
30 — 39 keV, which matches the first energy break ob- 



found that it significantly improves the Band-only fit by 
34 units of Cstat for two additional dof. Interestingly, the 



12 A BB spectrum peaks at an energy of about 3 times its tem- 
perature kT 
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TABLE 1 

Time-integrated spectral fit of GRB 120323A from To-0.016s to To+0.548s using the detectors n0, n1, n3 and b0 (see 

SECTION El. 









Standard 


Model 






Additional Model 


Cstat /dof 






B 


and, Compl 


or 2BPL 














E pcak 


a 
a 














kT or E 




B 


71 


-0.92 


-2.06 








- 


- 


- 


600 /470 




±5 


±0.07 


±0.02 
















2BPL 


40 
±2 


-1.22 
±0.03 


-1.96 
±0.02 


1024 
+ 282 
-213 




-5.35 
+ 1.75 








540/468 


B + Cut 


62 

±4 


-0.82 
±0.08 


-1.98 
±0.02 


764 
+ 308 
-450 


234 
+ 592 
-279 










551/468 


B + BB 


263 
+80 
-44 


-1.44 

+0.05 

-0.07 


< -2.36 












11.29 
+ 1.21 
-0.68 


568/468 


C+BB 


307 
+24 
-20 


-1.48 
±0.03 














11.72 
±0.50 


567/469 


B + C2 


239 
±25 


-1.45 
±0.03 


-2.49 

+0.15 

-0.28 








+2.45 
+0.57 
-0.51 




9.61 
+ 1.51 
-1.27 


568/467 


B + B2 


369 
+48 
-0.43 


-1.60 
±0.05 


< -10 








+ 1.64 
±1.05 


-2.46 

+0.19 

-0.29 


12.51 
+ 6.25 
-3.21 


549/466 


C + C2 


433 

+0.27 

-0.15 


-1.27 
±0.05 










+0.68 
+0.49 
-0.22 




38.93 
+ 10.3 
-11.3 


559/468 



tained with the 2BPL and B+Cut models. Further, the 
Epeak of the Band function is compatible with the sec- 
ond energy break at several hundred keV obtained with 
2BPL and B+Cut. 

However, with a B+BB model, we can only determine 
an upper limit for /3, which makes the Band function 
similar to a Compto nized function (a pow er law with an 
exponential cutoff - 
similar parameters 



Kaneko et al. 20061; both fits give 
peak and a. This is also evident 
from the fact that the combination of a Comptonized 
function with a BB (C+BB) which has only one dof dif- 
ference from Band, leads to the same Cstat improvement 
as B+BB. The Cstat change (ACstat) per dof between 
Band and C+BB is the largest of the tested models. The 
ACstat per dof differences between Band fits only and 
C+BB, 2BPL, B+Cut and B+BB are 33, 30, 25, and 
16 respectively. We discuss below the uniqueness of the 
various selected components in our fits. 

As shown in Figure [3] (bottom right panel) , the BB 
component appears like a hump in the low energy power 
law of the Band function. In order to better explore 
the intrinsic shape of this hump, we replaced the BB 
component with less constrained shapes such as another 
Band (B+B2) or a Comptonized (B+C2) function. With 
more parameters, the Band (i.e., B2) and Comptonized 
(i.e., C2) functions have no reason to mimic the shape of 
a thermal Planck function, which can be approximated 
with a Band function with a = 1 and a very steep (3, 
which is also equivalent to a Comptonized function with 
a power law index value of 1. In addition, due to repro- 
cessing of the thermal emission, a pure black body shape 
is very unlikely to be obtained. Interestingly, for both 
B+B2 and B+C2, the a values of B2 or C2 are positive 
and range between +0.60 and +3.00 and the spectral 
peaks of B2 and C are compatible with the temperature 
of the BB. This shows that the hump identified in the 
low energy power law of the Band function is compatible 
with a thermal origin and is adequatelly approximated 
with a BB component. 

Next, we tested the effect of the constrained curvature 
of the Band function by replacing the Band function with 



a smoothly broken power law (SBPL) with a free break 
scale, in the B, B+Cut and B+BB models. We obtained 
a slight improvement of 19, 10 and 10 units of Cstat, 
respectively, which does not impact drastically the fit 
results nor the model comparisons. Therefore, in the 
following we will use the Band function since with one dof 
less it is easier to derive well constrained fits, especially 
in shorter time bins and when several components are 
used. 

We note that the global shapes of the three favored 
models, 2BPL, B+Cut and B+BB (see Figure [3} are 
very similar and indicate that the time-integratea vF v 
spectrum is better fit with two bumps rather than the 
single one of the Band function. However, none of the 
tested models lead to a completely satisfactory fit of the 
time integrated spectrum. Whatever the model, system- 
atic patterns remain in the fit residuals, which are not 
distributed randomly around zero across the studied en- 
ergy range (see Figure [3]). This indicates that the mod- 
els tested are either not sufficient enough to describe the 
data, or that a possible strong spectral evolution during 
the burst leads to an unsatisfactory description of the 
spectrum with standard models when integrated over the 
whole burst duration. It is very likely that our models 
could present a good description of the data when in- 
tegrated over time scales encompassing periods with no 
or less spectral evolution. To investigate such possible 
issues as well as to follow the evolution of the spectrum 
during the burst, we performed time-resolved spectral 
analysis at shorter time scales, presented in Section [5] 

We compared the spectral results obtained with GBM 
to the Konus-WIND spectra to ensure that there is 
no major calibration problems. We analyzed GBM 
data from To+0.002s to Tp+0.256s, which is a simi- 
lar tim e interval to the one used in IGol cnctskii et al.l 
(2012) taking into account the propagation time be- 
tween the two spacecrafts (V. Pal'shin private commu- 
nication). U sing a power law w i th exp onential cutoff as 
proposed in Golenetskii et al.| ( 2012[ ), we obtained an 
E P oak value of 251^^ keV and a power law index value 
of -1.65+0.02. These results are compatible with those 
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reported in Golenetskii et al. (|2012|) (i.e., E pea k=331+5Q 



keV and index=-1.57±0.07). Such a crosscheck mini- 
mizes calibration issues between the two different instru- 
ments and reinforces our confidence in the goodness of 
the data set we used for this analysis. 

5. TIME-RESOLVED SPECTRAL ANALYSIS 

We perform time-resolved spectroscopy of 
GRB 120323A at two different time scales, and re- 



sults are presented in Section 5.1 and 5.2 Our goals are 
(i) to investigate the effects of possible spectral evolution 
during the burst, and (ii) to track the evolution of the 
various spectral parameters of the fit components, such 
as power law indices, break energies and temperatures. 

5.1. Coarse time-resolved spectral analysis 

First, we defined four broad time intervals based on 
the main structures identified between 20 and 150 keV 
in the light curves presented in Figure [T] The first 
time interval from To-0.018s to To+0.058s includes the 
first peak of the light curve, the second from To+0.058s 
to To+0.100s covers the most intense part of the sec- 
ond peak, and the time intervals from T +0.100s to 
T +0.174s and from T +0.174s to T +0.600s corre- 
spond to the decay phase of the light curve. We fitted 
each time interval with the models used in section [4] the 
fit results are shown in Table [H 

We find that while the 2BPL and B+Cut models 
showed evidence for a high energy cutoff in the time- 
integrated spectrum, they do not lead to the same con- 
clusion when fitting time-resolved spectra. The Cstat 
improvement obtained with these models compared to a 
Band function only is modest in all four time intervals. 
This leads to the conclusion that the cutoff measured in 
the time-integrated spectrum is an artifact due to the 
strong spectral evolution present in GRB 120323A. This 
result also demonstrates that the measured 'spectral cut- 
off' at high energies based on the extrapolation of GBM 
time-integrated sp ectral fits to the LAT ener gy range (as 
proposed in e.g., Ackermann et al. (2012)) should be 



interpreted with caution. 

B+BB or C+BB lead to the largest Cstat improvement 
compared to Band-only, with a ACstat between Band 
and C+BB of 46 and 18 units for 1 dof difference in 
the second and third time intervals, respectively. We 
note here that with Band-only fits, we find unusually low 
values of E pea k (a few tens of keV) in the first two time 
intervals with the highest intensity, while the third one 
corresponding to the global decay phase of the light curve 
has an E pea k around 400 keV. The low E poa k values of the 
first two intervals are accompanied with high (positive) a 
values. The second and third time intervals exhibit lower 
values for a with power law slopes steeper than -1.5. In 
contrast, with a two component fit, B+BB or C+BB, 
E pea k is shifted to higher energies. This is particularly 
obvious for the second time interval, for which E poa k is 
shifted from -40 keV to -500 keV. The addition of the 
BB to the Band function also leads to a lower value of 
a and (3. When the BB is replaced with a Band or a 
Comptonized function like in B+C2, C+B2 and C+C2, 
the parameters of this function are similar to those of 
the BB function with similar Cstat value. These results 
confirm the existence of a hump in the low energy part 
of the spectrum. 



5.2. Spectral analysis of hardness-ratio selected time 
intervals 

We showed in Section |4| and |5 . l| that spectral evolution 
within a burst can affectdramatically spectral fit results, 
potentially leading to a misinterpretation of the physics 
of the observed emission. 

In with section, we defined time intervals as short as 
possible to reduce effects due to spectral evolution but 
still large enough to be able to adequately fit at least a 
Band or a Comptonized function to the data. We pro- 
pose a novel method to determine the time interv als for 
time-resol ved spectroscopy. Similar to the idea by Scar- 
gle ( 1998 ) of characterizing flux variation with Bayesian 



statistics, which is referred as Bayesian Blocks method 
(BBM), we apply BBM to the evolution of the GRB light 
curve hardness ratio (HR), which is a good proxy to its 
spectral evolution. First, HRs were calculated between 
8 - 100 keV and 105 keV - 10 MeV for combined Nal 
and BGO data in a base bin, which is the finest possi- 
ble bin with at least 25 counts in each energy band to 
ensure the Gaussian statistics of HR. The errors in the 
HR were propagated from the count errors. Then, BBM 
was applied to the HR profile to find its change points. 
The prior for BBM is chosen in order to have a suffi- 
cient number of counts to perform spectral fitting and to 
avoid individual time intervals including too much spec- 
tral evolution. This analysis resulted in 12 time intervals 
which are used to generate the light curve presented in 
the right panel of Figure [T| We verified our ability to re- 
construct properly a Band function and a B+BB model 
in all these short time intervals with simulations follow- 
ing the procedure described in Appendix [A] 

We fitted each of t hese intervals with the models pre- 
sented in section 5.1 (i.e., B, C, B+BB and C+BB) even 
if the Cstat improvement were not statistically significant 
based on the number of dof differences for the various 
models. These results are presented in Table H] as well as 
in Figures [4] and [5] When one model was clearly not ade- 
quate to fit the spectrum in a time interval based on the 
Cstat value, it was not included in Table |4j When the 
high energy power law index of the Band function could 
not be constrained, we estimated its la upper limit. 

The Band model gives a satisfactory fit to the data in 
all the short time intervals. However, we can constrain 
the parameters of the B+BB (or C+BB) model in all 
these time intervals, and C+BB is statistically signifi- 
cantl^p*] better than Band alone in the three time inter- 
vals T +0.054s to T +0.080s, T +0.080s to T +0.094s 
and To+0.094s to Tq+0.132s with a Cstat improvement 
of 9, 25, and 51 units for 1 dof difference between mod- 
els, respectively. In principle this could be over-fitting 
our data, in which case we would expect the BB com- 
ponent to pick up random statistical fluctuations in the 
spectrum, with erratic changes of temperature and nor- 
malization from one time interval to the other. Instead, 
the temperature follows a constant cooling trend identi- 
cal to the evolution reported from the coarser time inter- 
val analysis. This observational result can be hardly ex- 
plained with random statistical fluctuations of the num- 
ber of counts in the various energy channels of the de- 
tectors. 

13 Based on the procedure described in Section |a| 
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Fig. 4. — Fine time-resolved v¥ v spectra using Band function only. The solid line corresponds to the model obtained with the best 
parameters from the fit, and the thin lines correspond to the 1 a uncertainty on the best fit. 
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Energy (keV) Energy (keV) Energy (keV) 

Fig. 5. — Fine time-resolved uF v spectra using B+BB. The blue lines correspond to the Band function, the red ones to the BB component 
and the black lines are the sum of the two components. The solid black lines correspond to the model obtained with the best parameters 
from the fits, and the thin black lines correspond to the 1 a uncertainty on the best fits. 




The three time intervals for which the addition of the 
BB leads to the greatest Cstat improvement correspond 
to the second peak of the light curve where the BB was 
also a statistically significant improvemen t in the coarse 
time-resolved spectral analysis (see Section 5.1 1. Figure[6] 
shows the vF v spectra resulting from the fit to the data 
during the time interval To+0.094s to To+0.132 s using a 
Band function (left) and B+BB (right). The systematic 
pattern observed in the residuals of the Band fit is clearly 
flattened when adding the BB. 

6. ONE OR MULTIPLE COMPONENTS ? 

In this section, we describe the spectral evolution 
within the burst resulting from the fine time resolved 
spectral analysis presented in Section [5~2} and discuss the 
best two spectral fit models, a single component (Band 
or Comptonized) and a two component scenario (B+BB 
or C+BB). Figure [7] shows the evolution of the parame- 
ters of the various spectral components with time. The 
blue curves correspond to the parameters obtained when 
fitting Band-only to the data, while the red ones are ob- 
tained when fitting B+BB. 

6.1. Single component: Band junction 

When fit with a Band function only, the E pca k of 
GRB 120323 A tracks the burst flux, especially the two 
peaks identified in the light curves above 20 k eV (see 
Figure 11 1) as also seen in previous burst spectra ( Ford et 
al. 1995). However, while E pea k is usually harder m the 
earliest high intensity peaks, especially in cases of sim- 
ply structured light curves like GRB 120323A, here the 
first peak is rather soft with values ranging between ^40 
and 200 keV, while the second peak is much harder, with 
Epeak values reaching ^600 keV. Each pulse exhibits an 
intrinsic hard to soft evolution. The evolution of a shows 
a striking discontinuity at ^To+0.094s. During the first 
intensity peak of the light curve, a is mostly positive 
with values above +0.2 in some cases, while during the 
second one, the a values drop below —1.5. Similarly, the 
values of /? are well constrained between —1.6 and —2.4 
until 0.094 s after trigger time, while only upper limits 
below —2.7 can be measured thereafter. We note that 
the discontinuity in the evolution of the parameter val- 
ues appears simultaneously for all the parameters of the 



Band function. Figure [4] shows the evolution of the Band 
function with time. 

Interestingly, Rmfit converges towards two different 
minima when fitting a Band function to the data in the 
time interval T +0.080s to T +0.094s (see Table [§. 



The lowest Cstat is obtained for a low E 



peak 



value and 



a high a one. The other local minimum is obtained for 
a much higher E pca k va lue and a very steep a. We will 
describe in section 6.3 the impact of this result when 



comparing the single component to the two components 



scenario. 



6.2. Two components: B+BB Model 

When a B+BB model is fitted to the data, the Ep 0a k 
of the Band function undertakes a global hard to soft 
evolution all across the burst with values evolving from 
>3 MeV during the beginning of the burst to ~30 keV 
during the late intensity decay phase. However, E pea k 
tracks strongly the light curve flux with an increase of the 
values from ^200 keV to ^800 keV corresponding to the 
flux increase phase of the second peak of the light curve. 
The values of a remain mostly constant around —1.3; 
only upper limits (always below —2.4) can be determined 
for p. 

The temperature kT of the BB component decreases 
linearly with time in log-log space from ~40 keV to few 
keV with a possible plateau or small reheating at the 
time of the second peak of the light curve. However, it 
is difficult to confirm this small feature since it could be 
simply due to a correlation between the parameters of 
the two components due to the fit process. 

Globally, the Band function and the BB component 
evolve independently. Figure [5] shows the evolution of the 
two spectral components with time. The reconstructed 
photon and energy light curves in the same energy bands 
as those used for the count light curves in Figure [TJ arc 
presented in Figures [8] and [9j respectively. The BB com- 
ponent contributes more than half the emitted energy 
between 20 and 150 keV during the first peak of the 
burst until ~T +0.080s. This contribution decreases 
very quickly to a few percent during the second peak 
of the burst. However, the BB remains an energetically 
subdominant component when computed from the time 
integrated spectrum over the entire GBM energy range 
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Fig. 7. — Evolution with time of the parameters of the spectral fits using the shortest time intervals. The blue lines correspond to the 
Band-only fit to the data and the red lines correspond to the B+BB fits. 



(8 keV to 40 MeV), where it only contributes about 10% 
of the total radiated energy. 

Finally, we replaced the BB with a Comptonized func- 
tion (C2) in the time interval with the highest statis- 
tical significance for the existence of the BB (i.e., from 
To+0.94s to To+0.132s), to investigate the shape of this 
low energy excess. C+C2 leads to the same Cstat values. 
The parameter kT = I5.6I3 5 keV of C2 is similar to the 
temperature kT = 8.4 ± 0.4 keV of the BB in this time 
interval and the index of the C2 function has a value of 
+0.1^0 4 which is similar with what is expected from a 
perfect Planck function (i.e., +1). As discussed in Sec- 
tion [HJ a thermal emission component with a low energy 
slope index around +0.4 is expected due to reprocessing 
of the photospheric emission, which is compatible with 
our data set. 



6.3. Comparison 

The most striking result when comparing Band-only 
fits with B+BB ones is the dramatic difference in the pa- 
rameters of the Band function. The strong discontinuity 
observed in the evolution of the spectral parameters of 
the Band function around Tq+0.094 s with Band-only fits 



does not exist when fitting B+BB. In the B+BB model, 
Epeak is systematically shifted towards higher energies, 
and both a and /3 are shifted towards lower values. In the 
Band-only scenario, a values exhibit large variations be- 
tween the first and second peak of the light curve. With 
B+BB model, a remains constant throughout the burst. 
Similarly, the f3 values vary during the burst with Band 
fits only, while the Band function high-energy power- 
law is always compatible with an exponential slope in 
the B+BB scenario. Therefore, in the B+BB scenario, 
the Band function can be replaced with a Comptonized 
function with no impact on the Cstat value of the fit; 
this replacement is only possible after ~To+0.094s in 
the Band-only fits. This explains the discontinuity in 
the evolution of the high energy power law indices, /3, of 
the Band function around T +0.094s when fitting Band- 
only to the data as presented in Figure [7| 

Fundamentally, a comparison between a Band-only to 
B+BB fits defaults to comparing different global spec- 
tral shapes. The former corresponds to a single peak 
spectrum in the v¥ u space while the latter results in a 
two-peak v¥ v spectrum with each peak evolving inde- 
pendently. In the B+BB scenario, the additional BB is 
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a significant component in terms of flux, especially be- 
tween 20 and 150 keV, where it contributes to more than 
half of the total emission. The statistical significance of 
the additional BB component depends, for instance, on 
the energy separation of the Band and BB i/F^ peaks as 
well as on the relative intensity of the two components. 

In the next two paragraphs, we report the possibility 
for the low energy hump, which is well described with 
the BB component in the B+BB scenario, to be driving 
the Band function shape when fitting a Band function 
alone to the data in the first peak of the light curve (i.e., 
<T +0.054 s). We then point out the strong similarities 
between the Band function shape from the Band-only 
fit with the BB and the Band function shape from the 
B+BB fit in the first and second peak of the light curve, 
respectively. 

During the first peak of the light curve (i.e., 
<To +0.054 s), the E poa k obtained from Band-only fits 
decreases constantly from —200 to —40 keV. Similarly, 
the BB in the B+BB fits cools constantly from —35 to 



—10 keV during the same period of time. Since the max- 
imum of the v¥ v spectrum of a BB with a temperature 
kT is — 3xkT, then the maximum of the v¥ v spectra 
of the BB resulting from the B+BB fits decrease from 
-105 to -30 keV. The evolution of the peak of the BB 
spectrum from the B+BB fits is then very similar to the 
evolution of the E poa k of the Band-only fit. The positive 
values of a resulting from Band-only fits during the first 
peak of the light curve are also similar to the positive 
low energy slope of a Planck function. With both E pea k 
and a, the Band function from Band-only fits and the 
BB from B+BB fits have very similar shapes for the low 
energy part. The difference between the Band function 
and the BB appears in the values of f3. While a Planck 
function has very steep high energy spectral slope, the /3 
values of the Band-only fits are high during the first light 
curve peak (> —2.4). During the first peak of the light 
curve where the BB would be the most intense, this com- 
ponent could be driving the spectral shape when fitting 
Band alone to the data. This is illustrated in Figure 10 
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Fig. 9. — Energy flux evolution when using B+BB model in the same energy bands as the count light curves presented in Figure [I] For 
each figure, the top panel corresponds to the total energy flux (i.e., B+BB), the middle panel correponds to the energy flux per component 
(i.e., Band in blue and BB in green) and the bottom panel corresponds to the contribution of each component to the total energy flux in the 
corresponding energy band (i.e., Band in blue and BB in green). Note that the \-a uncertainties may be over estimated since covariancc 
matrices resulting from the spectral fit process are not used. 
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(left panel) where Band and B+BB fits are overplotted 
in a time-interval included in the first peak of the light 
curve. 

During the second half of the burst (i.e., >To+0.094 s), 
the Band parameters, E poa k, a and (3 obtained either 
with Band-only or B+BB are very similar. The E pca k of 
Band-only fits decreases from ^500 to ^20 kcV during 
the second peak of the light curve, and the temperature 
of the BB from B+BB cools from ~10 keV down to ~4 
keV, which corresponds to a vFu peak decreasing from 
~30 to ~12 keV. Thus, conversely to what is observed 
during the first peak of the light curve, during the second 
peak, the Band function shape resulting from Band-only 
fits is very different from the shape of the BB resulting 
from B+BB model but consistent with the Band compo- 
nent evolution obtained with B+BB fits. This is again 

(right panel) where Band and 
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illustrated in Figure 
B+BB fits are this time overplotted in a time-interval 
included in the second peak of the light curve. 

In section [6+| we pointed out that Rmfit converges to- 
wards two different minima when fitting a Band function 
to the data in the time interval T o +0.080 s to T o +0.094 s 
(see Table |4j). The best fit is obtained with a low E pea k 
value and a high a one (i.e., Bl) ; the other local mini- 
mum has a much higher E poa k value and a very steep a 
(i.e., B2). It is interesting to notice that with their mea- 
sured spectral parameters, Bl and B2 mimic the Band 
and the BB components of the B+BB fit, respectively. 
Before To+0.080 s, the fit of a single Band function to the 
data would mimic the BB component of the B+BB fit, 
while after To+0.094 s, the Band-only fit would mimic 
the Band component of the B+BB fit. This supports 
the hypothesis that the BB component from the B+BB 
model is very intense during the first part of the burst 
and would be driving the Band-only fit, while when the 
intensity of this BB component decays during the second 
part of the burst, the Band-only fit would be drive n b y 



the Band function of the B+BB scenario (see Figure 10 ) 
6.4. Light curve peak overlapping scenario 

We cannot completely exclude that the two humps de- 
tected in the vYv s pectr a when using the coarse time 
intervals (see Section 5.1) are an artifact due to spectral 
evolution when the two peaks of the light curve over- 
lap. Lets consider the time interval from To+0.058s to 
To+0.100s from Table[3]where the BB is statistically the 
most significant. This time interval includes the decay 
phase of the first peak of the light curve and the intense 
part of the second one. When fitting C2+C or C+B2 
to the data (see Table [3]), the component with the low- 
est Ep ea k has spectral parameters (i.e., both E pca k and 
a) compatible with the trend reported in Section 5.2 
when fitting Band-only or Compt to the data up to 
T +0.094s (see Table E|. Similarly, the spectral param- 
eters (i.e., both Ep Ca k and a) of the component with 
the highes t E p ca k are compatible with the trend reported 
in section |5.2 when fitting Band-only to the data after 



T o +0. 094s (see Table [2]). In addition, the flux of the 
component with the lowest E pca k is compatible with the 
decaying flux of the first peak of the light curve when fit- 
ting Band-only or Compt, and the flux of the component 
with the highest E pea k is compatible with the peak in- 
tensity of the second peak of the light curve when fitting 
the simplest models. Thus, the spectrum in this interval 
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Fig. 10. — i/F„ spectra of GRB 120323A in two time intervals, 
T +0.022s to T +0.044s (top) and T +0.022s to T +0.044s (bot- 
tom), included in the first and the second peak of the light curve, 
respectively. The solid blue and red lines correspond to the best 
fits using the Band function only and the two humps model (i.e., 
B+BB), respectively. The dashed lines correspond to the 1-tr un- 
certainties on the best fit. In the B+BB scenario, the low energy 
hump correspond to the BB component, while the high energy 
one correspond to the Band function. In the time interval from 
To+0.022s to T +0.044s, the Band function of the Band-only fit 
mimic the shape of the BB component from the B+BB model 
(i.e., low energy hump), while in the time interval from To+0. 022s 
to To+0. 044s, the Band function of the Band-only fit mimic the 
shape of the Band component from the B+BB model (i.e., high 
energy hump). This explains the low values for Epeak as well as 
the discontinuities in the evolution of the spectral parameters of 
the Band-only fits observed in Figure [7] 



could be described as the superposition of the end of the 
first structure of the light curve and the beginning of the 
second one. 

However, this scenario cannot explain all the observa- 
tions reported in this article. For instance, the super- 
position of the two peaks of the light curve can hardly 
explain the possibility to fit two components with simi- 
lar intensity at the very beginning of the burst where the 
contribution of the second peak of the light curve should 
be very weak. Further we note that a two-component fit 
at the beginning of the burst cannot be due to random 
fluctuations because of the monotonic trend of the BB 
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temperature (see Tableland bottom right panel of Fig- 
ure uh. It is also difficult to explain the evolution of the 
spectral parameters obtained when using the simplest 
models such as the sharp discontinuity observed simul- 
taneously for all parameters or the resulting Luminosity- 
Epcak relation described in section [7] 

7. E^ST rL BAND RELATION 



Golenetskii et al. (1983) reported for the first time a 
hardness-intensity correlation during the prompt emis- 
sion of GRBs observed in th e Konus experiment on the 
Venera 13 and 14 spacecraft. Borgonovo & Ryde (2001) 
confirmed this correlation in a sample ot BATSE GHBs 
extending it to the evolution of the Band E pca k during 
the burst and its corresponding Luminosity. | Guiriec "et| 
al. (20101 showed that the correlation between the Band- 
runction Ep ea k evolution and the count light curve in- 
creased with the energy range used to define the light 
curve for three short and bri ght GRBs observed with 
GBM. Ghirlanda et al. (2011a I extended this result to a 
sample of 13 short GRBs detected with GBM showing 
the correlation between the Band function luminosities 
and their instantaneous E pea k values, and then a lso to 
a samp le long GBM GRBs in [Ghirlanda et al | pOlO ' 



2011b). More recently, Lu et al.J l 20121 ) reported a sum 
lar analysis on a large sample ot 62 bright GBM GRBs. 
In Figure 1 1 (top left panel) , we plot the energy flux of 



the Band function (i.e., FP and p] between 8 keV and 40 
MeV versus the Band v¥ v peatenergy, E pea k,i, of a sin- 
gle Band function fit (thin line) or a B+BB (thick line) 
to the time resolved spectra of GRB 120323 A (blue) and 
GRB 110721A (red). While no correlation is observed 
between FP and and E poa k.i when fitting Band-only to the 
data of GRB 120323A, a strong correlation emerges be- 
tween the flux of the Band function, FP and , and E pca k,i 
when fitting a model with two humps (i.e., B+BB). We 
notice the large errors on FP and when fitting only Band 
compared to B+BB. We would expect the opposite be- 
havior since with more parameters in the B+BB sce- 
nario we would also expect larger uncertainties, which 
should propagate accordingly to the errors of FP and . The 
E P eak,i-FP and relation obtained using the B+BB model 
is simi lar to the previously repo r ted results using Band 



alone J Borgon ovo fc Ryde 2001 Ghirlanda et al. 2010 
2011a|bTjL = u et al.||201 ! 2[ r~We~conclude that this beha 



ior also favors the existence of a hump in the low-energy 
power law of the Band function. 

The difference in the results between Band-only and 
B+BB fits is mainly the shift of Ep^^i to higher ener- 
gies, when using the latter model as well as the decrease 
of the Band function contribution to the total energy 
flux, since anot her c omponent is intense at low energies 
(see also section 6.2). Figures M and [9] clearly exhibit the 
strong presence (over 60% of the total flux) of the BB 
in the first pulse between 8—150 keV, while this con- 
tribution becomes 40% and less during the second pulse. 

(top left panel): 
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These results are reflected in Figure 
the Band fit data points compatible with the relation ob- 
tained with B+BB correspond to the second peak of the 
light curve, where the BB is statistically the most signif- 
icant but at the same time the least intense component, 

14 i being the index of a time-resolved spectrum within a burst. 



thus affecting the Band E pea k,i in the least. In contrast, 
the points of the Band-only fit that lie off the straight 
line relation correspond to the first peak of the light curve 
when the BB is statistically less significant but most in- 
tense and thus would affect E poa k4 the most. This again 
reinforces the two component scenario (B+BB) in the 
first peak of the light curve. Although less extreme, sim- 
ilar results are obtained with GRB 110721A (red lines) 
(see Figure 1 1 ) for which an inten se BB emission was also 
identified (A~xclsson et al. 2012). The relation between 
pBand &n( ^ E peaki i seems to be especially strong during 
the decay phase of individual pulses. Detailed analy- 
sis of several very intense bursts is necessary to assess 
if Epeakj tracks the energy flux during the rising phase 
of a pulse or if a hard to soft evolution of E pea k,i is re- 
ally observed during this phase. In the latter case, the 
E P eak,i-FP and relation would only exist during the decay 
phase of ind ivid ual pulses. 

(top right panel), we plot the E pea k,i- 
GRB 090227B 
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In Figure 

pBand re i ation f two short GRBs , 



and 090228, analyzed in Guiriec et ah] ( |2010[ ) (thin line) 
together with GRB 120323A (thick line^Tn these two 
short GRBs, a weak BB could also be present, but it 
has very little effect on both the measured Band func- 
tion flux as well as on E pea k. Band-only is then a good 
enough model for these two GRBs. The three GRBs lie 
along the same E pca k.i-FP and relation. Short GRBs are 
nearby events with a narrow redshift dis tribution: ^80% 
of them have z <1 ( Racusin et al.pbll ). Therefore, the 
Epeak, ;-FP and relation is expected to lead to a similar 
correlation between the luminosity of the Band function, 
LP and , and Ep ^ ; in the rest frame of the central engine 

suggesting a possible universal LP and -E p c cak \ relation for 
short GRBs. The average short GRB distances would 
then be consistent with the similar relation for the three 
bursts in Figure [TT] 

We now extend this analysis including several long 
GRBs. Figure [ill (center) shows the E^ r LP and re- 
lations for both (TRB 120323A and GRB 110721 A. We 
used a typical short GRB redshift of 1 for GRB 120323A. 
Berger (20111 reported tentative spectroscopic redshifts 
for GRB 110721 A at either 3.512 or 0.382 based on ab- 
sorption line observations with GMOS on the Gemini- 
South 8-m telescope; the former value is c onsistent with 
the po ssible redshift of ^3.2 reported by |Greiner et al.] 
( |2011 1. However, doubt remains on the real identifica- 
tionof the afterglow for GRB 110721 A. The data for 
the two GRBs line up perfectly showing a very strong 
correlation between E p e e s * kii and LP and only when using 
a redshift of ~3.2 for GRB 110721A. This suggests that 
the long GRB 110721 A is harder and more intense than 
GRB 120323A in the rest frame. The blue and red 
dashed lines correspond to the fits to E p ^ k r Lp and data 
with a power law for GRB 120323A and GRB 110721A, 
respectively. The best parameters of these fits with their 
1-er uncertainties are: 



t Band 
lj 120323A,i 



(1.57+ 1.26)10 50 (Ep^V;) 1 - 32 * 013 



erg s 



-,50 /^rcst U.34±0.19 g -l 



LfwVi - (1-65 + 2.64)10- (E^r 
The simultaneous fit of the two data sets with a power 
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Fig. 11. — Top left: Band function v¥ v spectral peak versus the flux of the Band function computed from 8 keV to 40 MeV (i.e., E pea k,i- 
F, relation) obtained when fitting Band (thin) and B+BB (thick) to the time-resolved data of GRB 120323A (blue) and GRB 110721A 
(red). Top right: Band function uF y spectral peak versus the flux of the Band function computed from 8 keV to 40 MeV (i.e., Epeai^i-Fi 
relation) obtained when fitting B+BB (thick blue) to the tim e-resolved data of GRB 120323A (blue) and Band to the time-resolved spectra 



of two other short GRBs analyzed in 



Guiriec et al. 



I 2010 ) (thin) . Center: E pea k,i-Fi relation in the rest frame (i.e., E; 



pcak.i 



-L; relation) 



obtained when fitting B+BB to the time-resolved data of GRB 120323A (blue) and GRB 110721A (red). The thick dash blue and red lines 
correspond to the fit using a power law to the data of GRB 120323A and GRB 110721A, respectively. The thin black line corresponds 
to the simultaneous fit of all data with a power law. Bottom left: Band function u¥ y spectral peak versus the flux of the Band function 
computed from 8 keV to 40 MeV (i.e., Epea^^-Fi relation) obtained over a sample of short and long GRBs. Bottom right: E pea k,i-Fi 
relation in the rest frame (i.e., E p ^ k ; -Li relation) obtained over a sample of short and long GRBs. The thin black line corresponds to the 
simultaneous fit of all data with a power law. 
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law corresponds to the solid black line given by: 



(1.60±0.65)10 au (Epeak^) 



M.33±0.06 erg g -l 



(1) 



The very good consistency between these three relations 
supports the possible universal behavior of the Ep C c s * k i - 
L Band rc i ation . 

In Figure 
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(bottom left), we added the data for 



long GRBs 05D 916C and 090926A from |Guiriec et a~ 



to the sample of the GRBs presented 
1'hese two additional GRBs have measured red- 



in pr eparation 
above 

shifts estim ated at 4.35±0.1 5 ( Greiner et al. 2009) and 
at ^2.1062 (Malesani et al. |2009p , respectively, and are 
fitted with a combination or three components, a Band 
function, a BB comp onent and an addition power law 
( |Guiriec et aL|2011b| ). We notice that the E peakii -Ff and 
relations tor trie various GRBs are shifted, with the high- 
est z GRBs being also the dimmest. Figure 11 (bottom 
right) shows the E p e s ^ k ; vs Lf and for the same sample of 
GRBs, assuming a redshift of 1.0 for all short GRBs. All 
data points now line up; a power law fit (solid black line) 
gives: 

jBand = ( Lg g ± 0.84)10 5 ° (E rCSt M 33±0.07 



Vaka; er § s ( 2 ) 

The color dashed lines correspond to the individual 
fit to the data of each burst with a power law, and the 
results are : 

Lo9 a o2 d 27 Ba = (1.26 ± 1.65)10 50 (Ep^) 1 - 29 "- 1 ^ s" 1 

t Band 
-^090228,1 ■ 

t Band 
^0809160,1 

t Band 
lj 090926A,i 



(1.28± 3.22)10 50 (E p ° c s a t k i ) 1 - 29±a36 erg B - 
= (1.50± 4.02)10 50 (E r p c c s a t k!i ) 1 - 41±a30 erg s" 



(1.50± 0.28)10 50 (E p c o s a t k i ) 1 - 34±0 03 erg s" 1 

The relationships described by equations [T and [2] are 
identical within errors and reinforce the possible univer- 
sality of the E p c c s ^ k r LP and relation across short and long 
GRBs. We note here that this GRB sample (Guiriec et 



al. 2010 in preparation! ) is limited; an analysis using a 
larger sample is the subject of a followup study. Our 
GRB data sample is also too limited to quantify dis- 
persion effects which could be due to multiple physical 
parameters like for instance the bulk Lorentz factor and 
the jet opening angle, or possible selection effects which 
could prevent the detection of possible outliers to this 
relation. 

re- 



pcak,i 



The presence of the E poak!i -Ff and and E 
lations in these six GRBs leads to multiple conclusions: 

• The E pcakj i-FP and relation seems to be intrin- 
sic across the time-resolved spectra of single 
bursts, and the E p c s ^ k r Lp and seems to be sim- 
ilar fro m burst to burst on a la rge sample of 
events (Ghirlanda et al. 2011a|b I although the 
slope slightly differs from what has been previ- 
ously reported. It can, therefore, not be only at- 
tributed to instrumental selection effects a s often 
suggested t o explain the s o called Amati (lAmati 
et al ||2002|) and Ghirlanda jGhirlanda et al.|)2004| ) 
relations (Ghirlanda et al l|2008| |Kocevski||2012j f, 
nor to effects of the redshift correction simultane- 
ously impacting LP and and E p ° s * ki , resulting in an 
artificial correlation. 



• It has been suggested that the correlation between 
the parameters of the Band function and its vF v 
peak energy is driven by the intrinsic correlation 
between the parameters of the Band functio n itself 
(|Massaro et al.|2008| [Goldstein et al.|2012b[ ). With 
GRB 120323A and GRB 110721 A, we have counter 
examples showing that fits to the data with Band- 
only do not lead to any correlation between E peakj i 
and Ff and , while B+BB fits do. If the correlation 
were mostly a model artifact, instead of physically 
driven, it should also be present when fitting Band- 
only to the data. 

• Contrary to the Amati-like relations, the E p c s * k r 

IjBand relation does not lead to a universal scenario 
for the central engine, but instead to the similar- 
ity of the relativistic jet evolution and radiation 
mechanisms dissipating the energy released by the 
central engine. It is in fa ct an extension of th e 
so called Yonetoku relation ( Yonetoku et al.|2004 ), 
which correlates for a sample of bursts, the peak 
flux of each burst (integrated over 1 s) with its 
corresponding E poak . This relationship is similar 
to the E p g S * k ; -LP and relation albeit with only one 
data point per GRB. 

• For GRB 120323A the E pcakii -Ff and relation holds 
only when B+BB is fit to the data. Any physi- 
cal interpretation must reproduce this relation to 
be viable, thus making the E p e e s * k ._Fp and relation 
a tool that discriminates between theoretical sce- 
narios trying to explain GRB prompt emission. 
In addition, GRBs deviating from this relation 
when fit with the Band function only, such as for 
GRB 120323A and GRB 110721A, may include ev- 
idence for a strong additional component such as 
BB. We note tha t the outlier GRBs from t he E peak - 
L relation from |Ghirlanda et al.] ( |2010[ ), namely 
GRB 080916C, GRB 090510 and GRB 090902B are 
GRBs known to exhibit strong addi tional spectral 
components to the Band fu nction ( Guiriec et al. 



2010 2011b in preparation). Our sample also in- 



clude GRB 0809I6C and we showed that with a de 
tailed spectral analysis, GRB 080916C is perfectly 
consistent with our new E p ® s * k ._Fp and relation. 

• Yonetoku, Amati and Ghirlanda relations exhibit 
dispersion effects which could eventually be re- 
duced using a similar multi-component spectral 
analysis as the one proposed in this paper. 

• Finally, since (redshift corrected) data of both 
short and long GRBs satisfy very similar E p c s * k r 

IjBand relations, a well-calibrated and dispersion 
corrected formula could eventually be used to es- 
timate redshifts for GRBs in the absence of multi- 
wavelength follow-up observations. Such estimates 
would only require a relatively intense GRB to en- 
able accurate time-resolved spectroscopy. 

Beyond the previous conclusions, if the universality of 
this new Ep^ r LP and relation is confirmed, then it will 
open new perspectives for the development of future in- 
struments. In many cases such as in the gravitational 
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wave research field, the redshift of a GRB is a crucial 
quantity to measure and requires a complex chain of 
operations consisting of repointing telescopes at various 
wavelengths to the source. However, the initial localiza- 
tion is often not good enough to initiate the process at 
all. Here, the E p c s * k r LP and relation would allow a red- 
shift determination only from the study of the spectral 
evolution in the gamma ray emission of GRBs. Thus a 
large GBM-likc instrument with high sensitivity would 
be ideal to determine GRB redshifts. 



8. INTERPRETATION 

Gamma-ray bursts are associated with ultra- 
rclativistic outflows ejected by a newborn compact 



source (see e.g., Piran 2004) and their prompt emission 



is very likely due to internal dissipation within the ejecta 



( |Sari fc Piran|1997| . Assuming that GRB 120323A were 
a, standard short GRB with an intense BB component, 
it is tempting to associate this component to the pho- 
tospheric emission produced by the relativistic outflow 
when it becomes transparent at large distances from 
the central engine. Without any additional dissipation 
process at the photosphere, the predicted photospheric 
spectrum is indeed close to a BB with two main 
modifications : (i) the low-energy slope is affected 
by the complex geometry of the photosphere, which 
leads to a photon slope a ~ +0.4 instead of a = +1 



for an exact Planck fun ction (Goodman 1986 Pe'er 



2008[ |Beloborodov| |2010[ ); (ii) the observed peak can 
be broadened it the temperature evolves on a timescale 
which is shorter than the time interval used for the 
spectral analysis. The second effect should be limited in 
GRB 120323A as its brightness allows a refined analysis 
with time bins shorter than the duration of the two 
main pulses. 

The spectral analysis presented above is not sensitive 
to the precise value of the low-energy spectral slope of 
the low energy hump (i.e., BB component). While data 
require a component with a positive low energy spectral 
slope compatible with a Planck function shape to de- 
scribe the low energy hump, a modified BB component 
with a low energy spectral index around +0.4 does not 
affect the results and is perfectly compatible with the 
data as well (see |6.2[) . Then, by assuming that the BB 
component in GRB 120323A is a thermal component of 
photospheric origin, it is possible to put some constraints 
on three important physical parameters: the radius Rq at 
the base of the flow, the Lo rentz factor T, and the photo 



spheric radius i? D h (s ee e.g., Daigne fc Mochkovitch 2002 



Pe'er et al. 2007 and HascoeFet al. 1 2012c) tor a gener 



alization to the case ot magnetized outflows). These are 
given by: 







/nt 


3/2 




X 







_DlR_ 

2(1 + z) 2 



a T (l + z) 2 D L F BB 1-R 
m D c 3 1Z R 



(3) 



x [(I + ^/ntP 1 ^) 



Rph 



fT'j/ 



DlF BB K 3 



1 - R 



16m p c 3 



(l + zf 



R 



x [(I + ^/ntP 1 / 4 ,(5) 



where z and -Dl are the redshift and the luminosity dis- 
tance of the source, Fbb is the measured flux of the BB 
component in a given time bin, R = F BB /F to t is the ratio 
of the flux of the BB component over the total flux, and 
1Z is computed from _Fbb and the measured temperature 
of the BB component by 



TZ = 



( Fm 

Ut, 4 



BB 



1/2 



(G) 



In addition to these quantities that can be directly 
measured, there are two unknown parameters related 
to the GRB physics, the ratio cx//nx and the product 
(1 + cr) /nt, where ex is the fraction of the initial energy 
released by the source which is in thermal form (the 
initial fraction of magnetic energy is 1 — ex), /nt is the 
efficiency of the dissipation mechanism responsible for 
the non-thermal component observed in the spectrum, 
and a is the magnetization of the relativistic outflow at 
the end of the acceleration process. 

Unfortunately, the redshift of GRB 120323A is not 
known. We assume z = 0.5, which is a typical value for 
a short GRB. The tendency is that a lower redshift will 
reduce the constraints derived below but we checked that 
our conclusions are unchanged for 2 = 0. 1 or 2 = 1. Us- 
ing the res ults of th e sp ectral analysis (B+BB) presented 
in sections |5.1| and |5.2| we measure F BB , R and 1Z and, 
using equations (|3fl5j), we obtain Rq , T and i? p h listed 
in Table [2] The emciency of the photospheric emission 
/x = Lph/E, where L p h is the luminosity of the photo- 
sphere and E the injected energy flux in the relativistic 
ejecta, can be compared to the efficiency of the dissipa- 
tive mechanism responsible for the non-thermal emission 
/nt = Lnt/E, where Lnt is the non-thermal luminosity. 
Using the formulae above, we find 



h_ 
/nt 



R 



1 - R 



0.01 - 0.1 



(7) 



Since the ratio i?/(l — R) is within the range 0.01-0.1, 
the non thermal dissipative process is dominant in the 
case of GRB 120323A. 

For a pure fireball, the acceleration is thermal, so that 
ex = 1 and a = 0. Then, if the dissipation mechanism 
responsible for the non-thermal component has a very 
high efficiency (/nx — 1), the values listed in Table [2] 
give direct estimates of Rq, T and R p h- They are m 
good agreement with typical values expected for GRBs, 
except for the initial radius Rq which seems too large. 

If the jet opening angle is 9- 3 and the size of the ini- 
tial region, where the outflow is launched, is £, then 
i?o — If short GRBs are associated with the merger 
of a NS+NS binary system, the expected central engine 
is an accreting highly rotating black hole with a mass 
of 2 — 3 Mq . Then, the radius of the innermost stable 
orbit, which can give an estimate of £, is of the order of 
8—13 km (we assume a ~ 0.8 for the black hole spin). 
This leads to R < 700 km if 0j > 1°. The value of R Q 
obtained in Table [2] is a factor 10 — 10 3 larger. Even 
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TABLE 2 

Values of the initial radius, Lorentz factor, photospheric radius and total (isotropic equivalent) injected power derived 
from the results of the time-integrated and time-dependent spectral analysis, using equations |3j|5} . 



Time interval 


Initial radius ( 


cm) Lorentz factor Photospheric 


radius (cm) 


Total power (erg/s) 


T s tart 


T s top 


-Ro x [/nt/ct]" 


rx[(i + a)/ NT f/ 4 R p] 


,x[(l + 


c) /nt] 1/4 


E x / NT 








Time-integrated spectral analy 


SIS 








-0.016 s 


0.548 s 


4.2 x 10 8 


162. 


3.4 x 10 12 


2.5 x 10 52 






Time- dependent spectral analysis: 


4 bins 








-0.018 s 


0.058 s 


1.1 x 10 9 


172. 


7.3 x 10 12 


6.3 x 10 52 


0.058 s 


0.100 s 


2.9 x 10 9 


176. 


1.0 x 10 13 


9.5 x 10 52 


0.100 s 


0.174 s 


7.5 x 10 8 


138. 


8.8 x 10 12 


4.0 x 10 52 


0.174 s 


0.600 s 


2.2 x 10 s 


85. 


5.0 x 10 12 


5.3 x 10 51 






Time- dependent spectral analysis: 


12 bins 








-0.018 s 


-0.006 s 


4.4 x 10 7 


330. 


7.0 x 


10 11 


4.3 x 


10 52 


-0.006 s 


0.006 s 


9.1 x 10 8 


258. 


2.2 x 


10 12 


6.6 x 


10 52 


0.006 s 


0.014 s 


3.7 x 10 9 


218. 


3.5 x 


10 12 


6.1 x 


10 52 


0.014 s 


0.022 s 


3.2 x 10 9 


202. 


3.9 x 


10 12 


5.5 x 


10 52 


0.022 s 


0.044 s 


1.2 x 10 10 


161. 


5.8 x 


10 12 


4.1 x 


10 52 


0.044 s 


0.054 s 


1.6 x 10 10 


136. 


8.3 x 


10 12 


3.6 x 


10 52 


0.054 s 


0.080 s 


7.2 x 10 9 


157. 


8.5 x 


10 12 


5.6 x 


10 s2 


0.080 s 


0.094 s 


3.5 x 10 9 


167. 


1.1 X 


10 13 


8.8 x 


10 52 


0.094 s 


0.132 s 


2.6 x 10 9 


145. 


1.1 X 


10 13 


5.8 x 


10 s2 


0.132 s 


0.178 s 


2.0 x 10 9 


98. 


1.3 x 


10 13 


2.0 x 


10 52 


0.178 s 


0.250 s 


2.0 x 10 8 


96. 


6.6 x 


10 12 


TO x 


1Q 52 


0.250 s 


0.700 s 


2.7 x 10 8 


76. 


3.7 x 


10 12 


2.8 x 


10 51 



if higher black hole masses can be expected for NS+BH 
mergers, an initial size Rq above 1000 km seems quite un- 
realistic for most theoretical models of short GRB central 
engines. However, the assumption /nt — 1 is quite ex- 
treme. The main mechanism to dissipate energy above 
the photosphere for an unmagnetized outflow is the ex- 
traction of kinetic energy by internal shocks. T his pro- 
cess is known to have a low efficiency /nt ;$ 0.1 (Daigne 

& Mochkovitch 1998). As Rq oc (/nt/ct) 3 ^ 2 , a realistic 
efficiency /nt 



0.05 also leads to a lower value of Rq (by 

10 -2 ) which is in much better agreement 

The im- 



a factor /^? 

with theoretical expectations for short GRBs. 

pact on the other quantities is weak (a factor / N t^ 4 — 2) ■ 
The Lorentz factor is found in the range 160-660, in good 
agreement again with the theoretical expectations, espe- 
cially the constraints obtained from the 77 opacity ar- 
gument ( for recent discussion s in light of Fermi - LAT re- 
sults, see IRacusin et al.||20"TT] |Zhao et aL"l|2011| |Hascoet| 



et~aTl|2012bp . 



It is interesting to compare the value of the pho- 
tospheric radius, which is found to be in the range 
1.5 x 10 12 -2.7 x 10 13 cm, with the radius where the in- 
ternal shock phase starts, which is given by 



2T 2 ct v 



0.01 s 



(8) 



where t V ar is the minimum variability timescale in 
the initial distribution of the Lorentz factor of the 
outflow. Assuming /nt = 0.05, it is found that for 
short timescale variability £ var = 0.01 s, the ratio 
i?is/i? p h is in the range 2-200 except for bin 10 where 
Ri s /Rph ~ 0.97. As most of the non-thermal flux seems 
to be associated with variability on timescales larger 



than 0.01 s, the observed values are therefore fully 
compatible with the scenario where the non-thermal 
emission is due to internal shocks above the photosphere. 

Most GRBs, however, are not compatible with the 
simplest scenario where ex = 1, as is the case for 
GRB 120323A. At least in long bursts, it seems that 
no thermal component can usually be detected. In GRB 
100724B, where an additional BB component was found 
in the spectrum in a similar way as in GRB 120323A, 



this a dditional component was weaker (Guiriec et al. 
2011a). To reconcile most bursts with the standard GKB 



scenario, it is necessary to assume that ct < 0-1 (see 



e.g.JDaigne fe Mochkovitch|[2002^ ,|Zhang fe Pe'er|[2009 
Guiriec et al.|2011a||Hascoet et al.|2012a|c[ ). In this case 



most of the energy initially released by the source is in 
magnetic form rather than thermal and th e jet acceler 
ation can be magnetically driv en (see e.g. 
E1|1994( |Daigne fc Drenkhahn||2002l |Vlaha' 



Begelman fc 
kis fc K onigl 

20031 IKomissarov et al ||2009| |Tchekhovskoy et al.||2010| 



Komissarov et al.||2010| IGranot et al.J|2011|). This leads 



ioy 

if: 



to two possible scenarios : 

• The magnetization of the outflow at large distances 
from the central engine is still large (er > 1). Then 
internal shocks cannot form and the best candidate 
for the mechanism responsible for the non t hermal 

emission is magnetic reconnection (see e.g. Spruit 

et al.|[200Tl |Lyutikov fc Blandlmdl|2003l [Zhang fc 
Yan||20TT| K 



• Most of the initial magnetic energy is converted 
into kinetic energy (efficient magnetic acceleration) 
and the magnetization at large distances is low 
((7 < 0.1). Then, as in the standard fireball model, 
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internal shocks are the best candidate for the non- 
thermal mechanism. 

Does GRB 120323 A favor one of these two possibilities? 
A possible indication is given by the condition that the 
initial radius Rq is expected to be small for short GRBs. 



As illustrated in Figure 12 the condition Rq < 100 km 



leads directly, from the values listed in Table j2j to low 
values of the non thermal efficiency, /nt ^ 0.008-0.3 
if ex — 1, and even lower values of /nt if £t is much 
lower. Such efficiencies are in the range expected for 
internal shocks but are rather lower than the usually 
considered range for magnetic reconnection. A high effi- 
ciency /nt = 0.5 for magnetic reconnection would lead to 
Rq 1000 km for most bins, and even Rq > 10 4 km for 
some bins. Such a large value of the initial radius seems 
quite challenging for most theoretical candidates of the 
central engine of short GRBs, especially NS+NS mergers. 
Therefore, the first scenario (large magnetization a > 1 
at the end of the acceleration of the flow, leading to mag- 
netic dissipation as the dominant process) seems the least 
probable for GRB 120323A: except for especially ineffi- 
cient acceleration mechanisms, where most of the initial 
thermal energy is converted into magnetic rather than 
kinetic energy, this scenario would require a low val ue of 
6t to get a large final magnetization (see Figure 12 1 and 
then imply a very inefficient dissipative process. Even in 
the second scenario (efficient magnetic acceleration lead- 
ing to a <C 1 at large distance), a low initial thermal 
fraction ex ^ 0.1 in GRB 120323A implies a really low 
efficiency for internal shocks, typically /nt ^ 0.01 (see 
Figure 12). GRB 120323A would then represent a case 
where 6t is a little larger than in most GRBs. For in- 
stance, £t — 0.5 implies /nt ^ 0.05. The best candi- 
date for the internal dissipation above the photosphere 
is therefore internal shocks, as the magnetization at large 
distance is expected to be low in this scenario. 

From this discussion, it appears the observations 
of GRB 120323A are compatible with the simplest 
GRB scenario where the relativistic ejecta are thermally 
accelerated and the non-thermal emission is produced by 
internal shocks. It is also compatible with the scenario 
where the fireball is initially magnetized and where 
most of the magnetic energy is converted into kinetic 
energy (efficient magnetic acceleration, a < 0.1) so 
that the dominant dissipative process remains internal 
shocks. GRB 120323A would represent a case where 
the initial magnetization is lower than in other GRBs 
(eT — 0.5-1 rather than 0.1 or lower). On the other 
hand, it is difficult to reconcile the data with a scenario 
where the outflow is magnetically dominated at large 
distance (cr > 1) and the non thermal emission is 
due to magnetic reconnection, unless this dissipative 
mechanism is much less efficient than usually considered. 

The spectral evolution observed in GRB 120323A can 
now be discussed in the framework of the pref erred sce- 
nario identified above. As shown in Section [5. 2| the spec- 
tral analysis based on the B+BB fits leads to a dramatic 
change of the low-energy slope a of the Band compo- 
nent, compared to the Band-only analysis. Instead of 
very steep values, often well above a = 0, the low-energy 
slope a is found for the B+BB analysis to remain for 
the whole burst in the range —1.9 — > —1.1. This is well 
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Fig. 12. — Constraints on the efficiency of the non-thermal 
dissipation process. Bottom: the maximum value of /nt given 
by the condition Ro < i?o,max = 100 km is plotted as a function of 
the initial thermal fraction ex during the energy release by the cen- 
tral engine, using the results listed in Tablep](red: time-integrated 
analysis; shaded area: time-resolved analysis ; the two limits cor- 
respond to bin 1 (-0.018 -> -0.006 s) and bin 6 (0.044-0.054 
s)). If the value of i?o,max is increased, all the lines are shifted 
towards higher values of /nt- A dashed horizontal line indicates 
the location of /nt — 0.1. Top: the magnetization cr at the end 
of the acceleration of the outflow is plotted as a function of et 
(thick solid line) , assuming no conversion of magnetic energy into 
kinetic energy, or of thermal energy into magnetic energy, and as- 
suming a perfect conversion of thermal energy into kinetic energy 
(i.e., a passive magnetic field frozen in the expanding fireball, with 
a final magnetization Cpassivc = (1 — £t)At)- The shaded area 
with <r < <Tp asa i ve corresponds to efficient magnetic acceleration 
models. The non-shaded area corresponds to scenarios where the 
magnetic field does not contribute to the acceleration and where, 
in addition, part of the thermal energy is converted into magnetic 
energy. Such scenarios are usually not considered in GRB mod- 
els. A dashed horizontal line indicates the limit a = 1 between 
magnetically dominated outflows (cr > 1, where the best candi- 
date for internal dissipation is magnetic reconnection, and matter 
dominated outflows (cr < 1), where internal shocks are the best 
candidate for dissipation. 



below the synchrotron slow cooling limit (a — —2/3), 
and well inside the predicted range for the fast cooling 
regime. The latter regime is expected during most of the 
prompt phase, i.e., —1.5 < a < —1, where a = —1.5 
corresponds to p ure fast cooling synchrotron radiation 
( Sari et aT]|1998[ ) and where steeper values are obtained 
it low-energy photons experience in verse Compton scat 



terings in the Klein Nishina regime j Dcrishev et al.|2001 
Bosnjak et al.|2009|[Nakar et ai.|2009||Daigne et al.|201l| 

Therefore, the spectral analysis of GRB 120323A based 
on the B+BB fits agrees well with the scenario where the 
non-thermal prompt soft gamma-ray emission is dom- 
inated by fast cooling synchrotron radiation of shock- 
accelerated electrons in internal shocks. A si milar re- 

f2011a j) in GRB 
RFT20323A of- 



sult was already found by Guiriec et al. 
100724B. Finding the same behavior in I 
fers a promising possibility to solve, at least partially, the 
so-called synchrotron death line problem (Preece et al. 



20 



Guiriec et al. 



1998] iGhisellini et al.|[2000| . 

It should also be noted that in the B+BB interpre- 
tation, the spectral evolution observed for the Band 
component follows the E poa k-Luminosity correlation 
observed in most GRBs (see Section [7]). This has been 
studied by several authors in the context of the internal 
shock model: the spectral evolution is governed by 
the dynamical timescale associated with the propaga- 
tion of shock waves which r eproduce successfully the 
E peak - Lumi n osity correlatio n) Daign e~ fc Mochkovitch] 
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Bosnjak et al . 2009 ; Asan o & Meszaros, 
Daigne et al.| |2011| |Bosnjak and Daigne| 
On the other hand, no similar strong correlation 
is found in GRB 120323A between the flux and the 
temperature for the BB component. As both quantities 
have a similar dependency with the Lorentz factor of 
the outflow, but a different dependency on E, this may 
indicate that not only the Lorentz factor, but also other 
parameters such as the total injected power, are variable 
during the relativistic ejection by the central engine. 
Finally, it should be noted that for the low efficiency 
/nt — 0.01 — 0.05 implied by the scenario discussed 
here, the isotropic equivalent total power of the outflow 
should be of the order of E ~ 6 x 10 52 -> 9 x 10 54 erg/s 
to reproduce the observed flux (from the values listed 
in Table |2j, leading to a total isotropic equivalent 
energy _E; so = (1.5 — 7.5) x 10 53 erg, which would favor 



a small opening angle 9\ 



1° for the true energy 



E ~ (Of/2) E iso to be consistent with the energy budget 
discussed for inst ance for the pop ular NS+NS merger 
scenario (see e.g., Aloy et al.|2005 1. Such a value of 0j is 
smaller than what is usually considered for short GRBs 
(iRuffert fc Janka||1999| jRosswog fc Ramirez- Ruiz||2002 
Rezzolla et al.||201l[ ). Obviously, this result is attectec 



by our choice of source redshift : the estimate of E^ 
reduced by a factor ~ 30 for z = 0.1. 

To illustrate the scenario where the prompt emission 
of GRB 120323A is associated with photospheric and 
internal shock emission in a (initially magnetized) 
relativistic outflow which is matter dominated at large 
distance (a <C 1), we have simulated a synthetic burst 
similar to GRB 120323A during the first - 0.2 s. We 
assume R = 300 km and ex = 0.5^ Other parameters 



are found in the caption of Figure 13 The photospheric 
emission is compute d using the method described in 
Hascoet et al. (2012c) and takes into account the modifi- 
cation of the spectral shape due to the complex geometry 
of the photosphere. The final photospheric spectrum 
is, therefore, more complicated than a BB but the 
value of the temperature and the flux are reproduced. 
The internal shock emission i s computed using the 
multi- shell model developed by |Daigne fc Mochkovitch| 
( fl998| and assuming for simplicity a pure fast-cooling 
synchrotron spectrum (i.e., with a = —1.5), without 
any correction for inverse Compton scattering in the 
Klein-Nishina regime. Therefore, the low-energy slope 
observed in the Band component of GRB 120323A is 
not perfectly reproduced. The resulting spectrum (t otal 
and separated components) is plotted in Figure |13[ 
The overall shape of the spectrum of GRB 120323A is 
reproduced. This is encouraging but this scenario should 
clearly be investigated in more details, especially to test 
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Fig. 13. — The internal shock + weak thermal photo- 
spheric emission scenario. This scenario is illustrated by a 
synthetic burst showing the same spectral shape as observed in the 
B+BB analysis of GRB 120323A. A redshift z = 0.5 is assumed. 
The outflow is ejected at i?o = 300 km with an initial thermal 
fraction et = 0.5. The duration of the ejection is 0.2/(1 + z) 
s. The acceleration is efficient so that a <g 1 at large distance. 
The initial Lorentz factor varies between 90 and 450 on a typical 
timescale 0.1/(1 + z) to produce two main episodes of emission. 
The isotropic equivalent kinetic power is assumed to be constant 
and equals 3.5 X 10 54 erg.s - 1 . The non-thermal emission from in- 
ternal shocks is computed assuming that a fraction e c = 1/3 of the 
dissipated energy is injected in relativistic electrons, so that the ef- 
ficiency is this case is /nt = 0.023. The observed time-integrated 
spectrum (over a duration 0.2 s) is plotted in black. The fast cool- 
ing synchrotron component due to internal shocks is plotted in blue 
and the thermal photospheric emission in red. As explained in the 
text, the spectrum of the latter is not an exact Planck function. 
The low-energy slope of the non-thermal component does not fully 
agree with the observations of GRB 120323A as we have assumed 
pure fast cooling synchrotron {a = —1.5) for simplicity (see text). 



if the observed spectral evolution of both components 
can also be reproduced. 

The discussion above has been limited to GRB mod- 
els which assume that the non thermal emission is due 
to internal dissipation in the relativistic outflow above 
the photosphere. There is however another theoretical 
possibility, where the whole spectrum would be of photo- 
spheric origin. The spectrum originating from the photo- 
sphere can be significantly modified if there is additional 



dissipation close to the p hotosphere ( Rees fc Meszaros] 
2005| |Pe'er fc Waxman||2005) . Such dissipation pro- 



cesses could tor instance be due to internal, or oblique, 
shocks, magnetic rec onnection (Giannios 20081, or col- 
lisional mechanisms (Beloborodov 2010). If the dissipa- 
tion produces a population of energetic leptons and a 
strong magnetic field, the original Planck spectrum can 
be modified by Comptonization, causing the spectrum to 
extend to higher energies, and by additional low-energy 
synchrotron photons, causing the spectrum to extend to 
lower energies. Depending on the conditions at the dis- 
sipation site, foremost the optical depth, the spectrum 
partly thermalizes again. These processes are capable 
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of produc ing a broad spectr um much resembling a Band 
function (Pe'er et al. 2006), but possibly showing steep 
low-energy slopes. The conditions that need to be met 
are that the energy given to the electrons should be com- 
parable to the energy in thermal photons and that a 
strong magnetic field exists . The details of the spectrum 



formation can be found in Pe'er et al. (20061; Vurm et 



al. (2011) and Ryde et al. (2011). This scenario seems a 
natural candidate to explain the results of the Band-only 
spectral analysis. As for the previous discussion, the ca- 
pacity of this scenario to reproduce the observed spectral 
evolution needs to be tested in details. In particular, one 
puzzling fact must be investigated furthermore: the tran- 
sition from the quasi-thermal (a^O) to the non-thermal 
(aj-1.5) spectrum at T o +0.094 in GRB 120323A. This 
may be a signature of a threshold for the dissipative pro- 
cess to occur at the photosphere su ch as suggested for 
instanc e in the collisional model by Bclobor odov et aT] 
(12011b). 



9. CONCLUSIONS 

We have presented here observational results and their 
associated theoretical interpretation of GRB 120323A, 
the most intense short GRB observed thus far with 
GBM. This GRB is especially bright below 150 keV. 
We associated the extreme intensity of the soft energy 
photons with the presence of a photospheric compo- 
nent. It is arguable whether the intensity of this event 
is entirely due to the existence of this component, de- 
tected for the first time due to the low energy range of 
GBM (starting at 8 keV) and possibly the vicinity of 
the source. GRB 120323A is either an unusually soft 
and intense short GRB or a regular short GRB but ex- 
hibiting an intense additional thermal-like component at 
low energy. Regardless of the origin of the component, 
GRB 120323,4 is a rare event. 

In summary, our observational analyses of the prompt 
7— ray emission of GRB 120323A have led to the follow- 
ing conclusions: 

• The presence of spectral evolution during a burst 
can create artificial features in the spectral shape. 
One should, therefore, be very cautious when in- 
terpreting time-integrated spectra. Time-resolved 
spectroscopy is required to remove these effects. 

• We determined that the spectra of 
GRB 120323 A (time-integrated and time- resolved) 
are better described with a double curvature 
spectral shape than with the single curvature of 
the Band function. The spectrum can thus be 
interpreted as consisting of two components, one 
being of thermal origin, compatible with a BB 
or similar shapes with steep slopes (i.e., >0) and 
the other produced by a non-thermal radiation 
mechanism. The former component is energetically 
subdominant compared to the non-thermal one. 

• The simultaneous fit of a thermal and a non- 
thermal component to the data dramatically 
changes the shape of the spectra of GRB 120323A. 
Using a single Band function we find that the 
time evolution of all spectral parameters exhibits a 
pronounced discontinuity, whereas all parameters 



evolve very smoothly in the two components sce- 
nario. In the latter scenario, the thermal compo- 
nent, whose shape is compatible with the expected 
shape of the photospheric emission of a relativistic 
jet, is most intense at the beginning of the prompt 
emission with a constant cooling trend thereafter, 
which closely follows its intensity decline. The pa- 
rameters of the non-thermal component are com- 
patible with fast cooling synchrotron emission. 

• Intriguingly, no correlation is found between the 
Epcaik i an< ^ ^ e luminosity of the burst, L;, when 
using a Band-only fit, while a strong correlation 
is obtained between L; and the Ep° s * k ; of the non- 
thermal component, when both thermal and non- 
thermal emission are fit simultaneously. In the lat- 
ter case, the GRB 120323A E^ s * k)i -Lj relation is 
perfectly consistent with those reported for a larger 
sample of GRBs. This result reinforces the two 
component scenario and supports the physical ori- 
gin of this relation as well as the possibility to use 
it as a discriminator for the prompt emission mod- 
els. This relation could also eventually be used as 
a possible redshift estimator for cosmology. From 
our (limited) sample, we estimate 

L p a „d = (1.59 ± 0.84). 10 50 (E™ f ! i t k i ) 1 - 33±0 07 erg.s- 1 . 

However, a more detailed analysis on a larger sam- 
ple of GRBs is required to estimate the dispersion 
of this relation. 

Our theoretical interpretation leads to the following 
conclusions: 

• The single component spectral analysis of 
GRB 120323 A clearly favors a photospheric origin, 
due to the steep low-energy slopes that are not 
compatible with the synchrotron radiation from 
shock-accelerated electrons. Additional dissipative 
processes at the photosphere are however necessary 
to reproduced the shape of the observed spec- 
trum. It remains to be tested if this dissipative 
photospheric emission scenario can reproduce 
the observed spectral evolution and especially 
the transition from a quasi-thermal (o^O) to a 
non-thermal (aj-1.5) spectrum at T +0.094s. 

• The two component spectral analysis of 
GRB 120323A is compatible with the standard 
fireball scenario where the thermal component 
is associated with the photospheric emission of 
a relativistic jet that has been thermally accel- 
erated, and the non-thermal emission is due to 
synchrotron radiation from accelerated electrons 
in internal shocks above the photosphere. It is also 
compatible with an alternative scenario suggested 
by other bursts, such as GRB 100724B, where 
the initial magnetization of the ejecta is large and 
where most of the magnetic energy is converted 
into kinetic energy below the photosphere (effi- 
cient magnetic acceleration). GRB 120323A would 
however correspond to a case where the initial 
magnetization is lower and the thermal energy 
content larger than in previous cases, such as 
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GRB 100724B. Indeed, our analysis shows that 
a large initial magnetization would lead to a 
very low efficiency of the non thermal emission 
and therefore to an energy crisis for this burst. 
On the other hand, it is difficult to reconcile 
GRB 100724B with a scenario where the outflow 
is still highly magnetized at large radius and 
where the non-thermal emission is due to magnetic 
reconnection above the photosphere. This would 
lead to an injection radius where the outflow is 
launched much too large for most models of the 
central engine of short GRBs. 



• The physical interpretation of the two-component 
spectral analysis suggests that the same processes 
(photospheric emission + synchrotron radiation 
from electrons accelerated in internal shocks) may 
be at work in both short and long GRBs, but that 
the composition of the relativistic jet may slightly 
differ in the two classes, due to different progeni- 
tors. GRB 120323A data suggest an initial magnc- 
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tization which is lower in short GRBs. The detec- 
tion of a two component spectrum in other short 
GRBs would be necessary to confirm this result, 
which can shed light on the nature of the central 
engine and the jet acceleration mechanism in the 
different classes of GRBs. 
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APPENDIX 
SIMULATIONS 

We used Monte Carlo simulations to (i) validate our ability to reconstruct properly the spectra in our shortest time 
intervals where the statistical fluctuations could be critical, and (ii) determine the probability of one model being 
better than another. 



Technique 

For each simulated set (i.e., one time interval and one input model), we generated 30,000 synthetic spectra for each 
relevant detector, with each spectrum covering the same duration as the real source time interval. For each synthetic 
spectrum we sum in each energy channel the number of background counts estimated from the real data (by fitting 
a polynomial function to off-source time intervals and by extrapolating it during the source active period), and the 
number of counts expected from the theoretical input photon model (i.e., Band or C+BB with the parameters obtained 
when fitting the real data in this time interval) when folded through the detector response matrix used for fitting the 
real data. Poisson fluctuations are then applied in each energy channel to the sum of the signal and background counts. 
For each data set background is also simulated based on the real background fit. 

During the simultaneous fit of the data from all the detectors, a new background is simulated for each synthetic spec- 
trum and each detector by adding Poisson fluctuations in each energy channel to the previously estimated background. 
The rest of the fit process is the same as with the real data fits. 

We performed all simulations using the same version of Rmfit as the one used to fit the real data. 

Model comparison 

We computed the significance of the improvement of a C+BB fit over a Band-only fit in the time intervals To-l-0.080 s 
to To+0.094s and To+0.094s to To+0.132s, where the additional BB component improves the Band-only fit by 25 
and 51 units of Cstat, respectively. We defined a set of simulations using Band as the input model (i.e., null hypothesis) 
with parameter values being the central values obtained by fitting the real data. We then fitted the resulting synthetic 
spectra with both Band and C+BB. 

The resulting distributions of the Band function parameters obtained from the fits of the synthetic spectra have 
symmetrical and peaked shapes for E pca k and /3, while the distribution of their amplitudes is wider. In the cases of 
positive values of a in input, the distribution of a is less peaked, but still with positive values. Taking into account 
the errors in the parameters, the input parameters are adequately reconstructed within 2a. This indicates our ability 
to reconstruct a Band function in this time interval, given the parameters of the Band function, the level of the 
background and the signal strength. 

We then fitted the synthetic spectra with C+BB ; most of the 30,000 fits gave a higher Cstat value than those 
using Band. Although C+BB has 1 dof more than Band, C+BB is usually worst than Band alone with Cstat values 
larger than 100 units, besides a few cases were C+BB improves the Band-only fit by only a few units of Cstat. For 
instance, for the time interval To+0.080s to To+0.094s, the maximal value for ACstat computed between Band and 
C+BB is ^3 units from the simulation. We observe a ACstat of 25 units in the real data set, which indicates that the 
probability that C+BB is better than Band due to statistical fluctuations in the real data is very likely to be much 
lower than 3.3xKT 5 . 

We performed a similar analysis chosing C+BB as the null hypothesis. When fitted with C+BB, the parameter 
distributions are narrow and symmetrical showing the very good quality of the reconstruction when using this model 
for the two time intervals. When fit with a Band function alone, the resulting parameter distributions are perfectly 
compatible with what is obtained when fitting the real data with a Band function. For instance, for time interval 
To+0.080s to To+0.094s, E pea k is reconstructed around 50 keV with a dispersion of a few keV, a has values around 
and the distribution of /3 values is compact around -2. In addition, the distribution of ACstat measured between Band 
and C+BB from the synthetic spectra is narrow and peaks around the observed value (i.e., 25), which was not the 
case when choosing Band as a null hypothesis. This indicates that if C+BB is the correct description for the spectral 
shape of the real data, then we expect from the simulation that a fit of the real data with Band-only would lead to 
the observed results. Conversely, if Band is the real spectral shape, then based on the simulations we would expect 
different values for the spectral parameters when fitting C+BB to the real data than those measured. 

In both time intervals, C+BB is statistically significantly better than Band alone. 



Reliability of the reconstruction 

Using the same method as in Section 



A. 2 we verified our ability to reconstruct the observed spectra in two time 
intervals, T +0.022s to T +0.044s andTo+0.054s to T +0.080s, where the models Band and C+BB or C and 
C+BB, respectively, give similar Cstat values, but for which the spectral shape is dramatically different. In both 
time intervals, whatever is the input model, the input parameters are adequately recovered. We then confirm our 
ability to reconstruct the spectrum in such fine time intervals. We also confirm that the models are indistinguishable, 
and that when we fit the synthetic spectra with the model which is not the input one, we recover the parameters 
as obtained from the real data. This confirms that the two models are both possible options to describe the spectra 
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of GRB 120323A. From a statistical point of view, the model with the lowest number of dofs would be preferred. 
However, si nce t he more complex model makes more physical sense and is significantly better in some time intervals 
(see Section A.2|, both options are considered here. 
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TABLE 3 

Coarse time-resolved spectral fits of GRB 120323A using the detectors nO, n1. n3 and bO (see section 



Time 


Models 






Standard 


Model 






Additional 


Model 


Cstat/dof 






Band, Compt or 2BPL 


BB 


Compt 




r Band 




Tstart 


TBt ° P 


Parameters 


E peak 


a 


p 


E b 


E f 




C 






kT or B 




-0.018 


+0.058 


B 


73 


+ 0.10 


-2.16 








- 


- 




- 


506/470 








±5 


±0.15 


±0.03 




- 




















43 
±3 


-0.51 
±0.10 


-2.03 
±0.04 


810 
+ 167 

-320 




-2.83 
+ 0.41 
-61 


- 


- 




- 


514/468 








69 

±4 


+ 0.19 
+ 0.17 

-0.15 


-2.10 
±0.04 


718 
+ 1000 
-300 


700 
fix 


fix 


- 


- 




- 


504/468 








118 

±14 


-0.14 
+ 0.34 
— 0.23 


-2.23 
±0.04 


- 


- 


- 


- 


- 




10.50 
+ 1.11 
-1.07 


495/468 








105 
±11 


-0.40 
±0.12 


-2.24 
±0.04 


- 


- 


- 


+ 7.35 
+ 3.5 
-4.5 


- 




4.11 

+2.49 

-1.83 


492/467 


+0.058 


+0.100 


B 


46 


+ 0.07 


-2.01 
















558/470 








±3 


±0.21 


±0.02 




- 


















2BPL 


32 
±2 


-0.87 
±0.09 


-1.94 
±0.03 


1040 
+274 
-194 




-3.94 
+ 0.85 
-1.91 


- 


- 




- 


534/468 






B + Cut 


44 

±2 


+ 0.10 
+ 0.27 
-0.19 


-1.94 
±0.03 


836 
+ 632 

-252 


393 
+ 704 
— 269 


- 


- 


- 




- 


532/468 






C + BB 


556 
+ 79 
-62 


-1.35 
±0.05 


- 


- 


- 


- 


- 


- 




10.94 
+0.40 


511/469 






B+BB 


547 
+ 78 
-62 


-1.35 
±0.05 


< -2.35 


- 


- 


- 


- 


- 




10.93 
+0.40 
-0.38 


512/469 






C + C2 


570 
+ 64 
-51 


-0.85 
+ 6.38 
— 0.30 


- 


- 


- 


- 


0.02 
+ 0.45 
— 0.24 


- 




22.62 
+4.72 
-4.96 


506/468 








+ 70 
-44 


+ 6.39 
-0.27 


- 


- 


- 


- 


0.05 
+ 0.42 
-0.28 


< -10 




22.12 
+5.24 
-4.50 


506/467 


+ 0.100 


+0.174 


C 


405 
+ 57 
-45 


-1.55 
±0.03 


- 


- 


- 


- 


- 


- 




- 


457/471 






B 


401 
+ 57 
-44 


-1.55 
±0.03 


< -3.40 


- 


- 


- 


- 


- 




- 


458/470 






2BPL 


19.59 
+ 6.27 
-2.90 


-1.26 
+ 0.23 
— 0.18 


-1.67 
±0.02 


505 
+ 100 
-71 


- 
- 


-3.23 
+ 0.37 
-0.70 










447/468 






B+Cut 


38 


-0.87 


-1.68 


270 


400 




- 


- 




- 


447/474 






C + BB 


395 
+ 48 
-38 


-1.37 
±0.07 




- 


- 


- 
- 


- 


- 




6.67 

+0.59 

-0.51 


- 

440/469 






B+BB 


391 
+ 47 

-38 


-1.37 
±0.07 


< -3.5 








- 


- 




6.65 

+0.58 

-0.51 


441/469 






C + C2 


373 
+ 41 
-31 


-1.18 

+0.20 

-0.16 










-0.07 
+ 0.94 
-0.58 






13.06 
+6.16 
-4.30 


439/468 






C + B2 


378 
+ 43 
-35 


-1.24 

+0.70 

-0.10 










0.27 
+ 0.60 
-1.00 


-5.92 




11.13 
+ 18.90 
-2.59 


439/467 


+ 0.174 


+0.600 


C 


61 


-1.77 


















516/471 








±5 


±0.04 
























B 


62 


-1.79 


-4.42 
















516/470 








±5 


±0.04 
























2BPL 


13 

+ 1.14 
-1.81 


-0.87 
+ 2.92 
— 0.43 


-1.95 
±0.03 


127 
+ 31 
-24 




-2.73 
+ 0.20 
-0.34 










512/468 






B+BB 


72 
±11 


-1.68 
±0.10 


< -2.66 














3.463 
+ 1.60 
-1.06 


514/468 
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TABLE 4 

Fine time-resolved spectral fits of GRB 120323A using the detectors nO. n1. n3 and bO (see section [5^2 



Standard Model 



Additional Model 



Cstat/dof 



nd, Compt i 



Tgtart | Tgtop 



BB, Compt 
index I kT or Eq 



-0.018 -0.006 



P + BB 
C+BB 
B + BB 



5S3 _125 
249 + 90 

6277+ 9010 

.,,-+10325 

£ " ! ' ! '-4225 



123 
1760 



+28 



+ 794 
604 



, „ 7 +0.04 

n S6+ ' 21 
-°- S6 -0 .37 

47+ ' 43 
u *'-0.27 

-1.38 + 0.08 
-1.21 + 0.12 
-1.19~* 



., + 0.19 
-0.11 



+0 31+ 44 
+o.ji_ 34 

-1.12 + 0.11 



+ 0.14 
3 -0.18 



O + 0.09 
-0.13 



40.T.+ "-™ 
.+9.19 
-7.49 



40.47^ 



n + 2.63 
-2.43 



386/472 
372/471 
364/470 
365/470 
362/469 
362/468 

417/470 
422/469 



C + BB 
B + BB 



109 + 13 

680+ 401 
bs0 -198 
, + 317 
-225 



162 ~ r 



+0 17+ ' 30 
+"■"-0.33 

-1.09 + 0.15 

-0 96+ ' 44 
uaD -0.19 



24 14+ 1 - 92 
. + 2.22 
' -3.39 



28.2JV 



420/470 
426/469 
423/469 



C + BB 
B + BB 



C + BB 
B + BB 



98 + 12 

497+ 203 
* M '-120 
+ 218 
■144 



50 l n 



62 + 5 
330 



+69 
48 



-0.29 + 0.23 
-1.31+0.11 
-1.32 + 0.13 

+ 54 + - 28 
+ u -° -0.24 

-1.16 + 0.11 

12 + - 66 
U ' -0.38 



-2.28 + 0.11 
< -3 

-2.34+0.06 



, , r + 0.09 
" 30 -0.09 



■?n 11+1.88 
2U11 -1.72 
1 + 1.89 
-1.78 



20:2 r 



15.53 + 0.73 
11.77+1.50 



389/470 
385/469 
386/469 

475/470 
475/469 
473/468 



+ 0.080 +0.094 



+ 0.094 +0.132 



C+BB 
B + BB 



C + BB 
B + BB 



C + BB 
B + BB 
C+C2 



C + BB 
B + BB 
C+C2 



C + BB 
B + BB 



44+5 
+ 49 
36 
+ 56 
40 



231 
111 



45 + 6 

476 + 146 
'"-105 

499 + 149 

735 + 217 
73 -146 

722 + 211 
' 22 -150 

743+?!! 

«*t85 

45S + 66 
°-52 

«<tS8 

429 + 37 

237 ig 

236 + | 2 

214+ 32 
zl -23 

194+ 30 



, „ 7 ,+0.56 

+°' 75 _0.49 

1 12 + ' 22 
1 - i -0.17 

-0.24 + 0.07 

+0.29+° 2 | 

1 36 + - 08 " 
iJO _0.07 

-1.36 + 0.08 

-0.08 + 0-04 
-1.60 + 0.04 
-1.60 + 0.04 
1 36+ ' 10 
1 35 + ' 10 



-1.36 



. 08 
+ 0-55 
0.15 



-1.50 + 0.03 
-1.50 + 0.03 
-1.14 + 0.10 
-1.14 + 0.10 
-0.73"^ 



, + 0.27 
-0.23 



-1.62 + 0.05 

-1.62 + 0-05 

-1.34 + 0.15 

1 24+ ' 20 
" 1 " ! *_0.17 



-2.29 + 0.07 

< -2.40 
-2.12 + 0.03 

< -3.3 

-1.93 + 0.03 

< -2.53 



„ + 0. 45 
-0.24 



+ 0.48 
-0.37 



+ 0.33 
5 -0.93 



10.99 n 



. + 0.83 

-0.75 

9 07+ 1 6B 
BU '-2.06 



11.79 
11.81 



10.16^ 



7.89^ 



15.58 n 



4.126 
4.17 



+ 0.54 
0.50 

+ 0.55 
0.52 



. + 0.63 
' -0.56 

10 14+ ' 61 
lu - 14 -0.56 

11 34+ 1 0-9 
11,M -4 .76 



. + 0.40 
-0.37 

7 89 + 40 
''° -0.37 
, + 4.55 
-3.47 



+ 0.58 
0.71 
+ 0.50 
0.60 



354/470 
352/469 
353/468 

480/470 
471/469 
472/469 

471/470 
508/470 
508/471 
446/469 
446/468 
446/468 

437/471 
438/470 
387/469 
388/470 
383/468 

461/471 
461/470 
456/469 
454/468 



C + BB 
B + BB 



122 + i 4 
122 + 1 4 

124+12 
124+" 



-1.50 + 0.06 

-1.50 + 0.06 

-1.42 + 0.17 

-1 41 + ' 16 
14 -0.18 



3 59 + 2 14 
3sa -1.15 

3.59 + 1'?? 



372/471 
372/470 
371/469 
371/468 



+ 0.250 +0.700 



18 + 10 

19 + 9 
35 + 17 



-1 


f>4±0 


06 








-1 


92±0 


07 


< -2.72 






-1 


85±0 


14 






3 



y -0.73 



564/471 
564/470 
563/469 



